Quantum phase transition in the Fulde-Ferrell-Larkin-Ovchinnikov states of a two dimensional 
d-wave superconductor upon increasing the exchange field 
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We studied the phase diagram for a two-dimensional d-wave superconducting system under an in-plane mag- 
netic field or an exchange field. According to the spatial configuration of the order parameter, we show that there 
exists quantum phase transitions in which the uniform phase transforms to the one-dimensional Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state, and then to two-dimensional FFLO state upon increasing the exchange field. 
The local density of states are calculated and suggested to be signatures to distinguish these phases. 



The Fulde-Ferrell- Larkin-Ovchinnikov (FFLO) state was 
predicted several decades ago by Fulde and Ferrell QJJ] , and 
Larkin and Ovchinnikov ||2|] for the superconductor in a strong 
magnetic field, where the superconducting (SC) order param- 
eter varies periodically in space. While the occurance of the 
FFLO state requires very stringent conditions on the SC mate- 
rials, namely, the Pauli paramagnetism effect should dominate 
over the orbital effect [3], and the material needs to be very 
clean [4[]. As a result, this long thought of inhomogeneous SC 
state has never been observed in conventional superconduc- 
tors. 

For layered systems with an exchange field or a mag- 
netic field parallel to the SC plane, the orbital effect 
will be suppressed strongly due to the low dimensional- 
ity. Thus they could be strong candidates to look for 
the FFLO state. Actually, in the past decade, indications 
for possible FFLO state have been reported in the heavy 
fermion materials CeCoIns |5l la], organic superconductors 
/t-(BETS) 2 GaCl 4 Jl,M BET S)2FeCi4 H § and k-(BEDT- 
TTF) 2 Cu(NCS) 2 ©El- All of them are quasi two dimen- 
sional (2D) layered compounds. The experimental develop- 
ments have attracted renewed interest on the property of the 
FFLO state. Theoretically, the existence and the character 
of the FFLO state can be investigated through analyzing the 
Ginzburg-Landau (GL) free-energy function, which is valid 
at temperatures not too below the superconducting transition 
temperature. Another effective method is the Bogoliubov- 
de-Gennes (BdG) technique and it has been proven to be 
a powerful tool to study the inhomogeneous state and the 
local density of states (LDOS) self-consistently in the low- 
dimensional system. In fact, in the past, the FFLO state 
has been studied intensivel y b ased on the above two tech- 
niques d H ElJllM H. For a 2D system, it is 
somewhat established I14l ll5Tl that the order parameter has a 
2D checkerboard pattern for a superconductor with c/-wave 
pairing symmetry, and a ID stripe-like pattern for s-wave pair- 
ing symmetry 11151 Il6ll . On the other hand, it has been shown 
that in the presence of dilute impurities, the pattern of the 
FFLO state becomes ID stripe-like in a 2D a!- wave supercon- 
ductor ddl]. This implies that the 2D and ID FFLO states 
may be present in the system in different parameter region. It 
was also argued without a calculation that the H—T phase dia- 
gram for the 2D isotropic systems, regardless it is s or c/-wave 
pairing, should include both ID FFLO state and 2D FFLO 
states with square, triangular and hexagonal patterns IU2I1 as 



the exchange interaction becomes stronger. This result is sig- 
nificantly different from our previous understanding about the 
pattern of the FFLO state in 2D systems. Therefore it is of in- 
terest and timely to reexamine the formation of the FFLO state 
in a 2D superconducting system as a function of the exchange 
interaction. 

In the present work, we calculate the spatially distributed 
order parameter self-consistently based on the BdG equations. 
The whole H -T phase diagram is constructed. We verify nu- 
merically that the pattern of order parameters in the FFLO 
state for two-dimensional cZ-wave square lattice samples in 
presence of an exchange field is not always simply 2D. At 
zero temperature, the pattern changes from uniform to ID 
FFLO state, and then to 2D FFLO state as the strength of the 
exchange field increases. The periodicity of the ID and 2D 
FFLO states decreases as the exchange field increases. At fi- 
nite temperature, the ID FFLO state will transit to the uniform 
phase upon increasing the temperature. Thus near the SC tran- 
sition temperature, only uniform phase and 2D FFLO state are 
observed. The LDOS for the above states are also calculated 
and they provide definitive signatures for the above mentioned 
FFLO states. In addition, our numerical study indicates that 
for a 2D s-wave superconductor, the FFLO state is always ID 
like and no 2D pattern could be obtained, this conclusion is 
consistent with that of Ref. lfl5l [Trill and the result will not be 
presented here. 

We start from a phenomenological BCS-type model with 
the Zeeman splitting effect caused by an exchange field or 
in-plane magnetic field. On a two-dimensional square lat- 
tice with a pairing interaction V between the nearest-neighbor 
sites, the mean-field Hamiltonian leading to the d wave super- 
conductivity can be written as, 



+ I li j(A ij cJc f .+h.c.), 



(D 



where f (/ are the hopping constants and pt is the chemical 
potential, crh is the Zeeman energy term, caused by the 
interaction between the magnetic field and the spins, with 
cr = +1 representing for spin-up and spin-down electrons, re- 
spectively. The af-wave SC order parameter has the following 
definition: Ay = V^c/jC^ - .c,-j,C/f)/2. 

This Hamiltonian can be diagonalized by solving the BdG 
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FIG. 1: H-T phase diagram of the two-dimensional d-wave super- 
conductor in the parallel magnetic field. 



equations, 



y/Hy Ay V U > \_ / uM 



where ffy is expressed by, 



(2) 



(3) 



The SC order parameter and the local electron density «, 
satisfy the following self-consistent conditions, 

A <v = T £ ( « + « ) tanh ^)< W 



^ l< T l 2 /(£„) + 2 K/f 1 " ( 5 ) 




FIG. 2: (Color online) Plots of the order parameter A as a function 
of position for various Zeeman fields h. 



Here fix) is the Fermi distribution function. We define the fer- 
romagnetic (FM) spin order m, and the on-site order parameter 

Aj as, nij = n n - «, u ; A,- = l/4(Ay +jE + Ay_i - A w+J - A,-, H .). 
The LDOS is expressed by, 

pA.cS) = J][|< T | 2 <5(£„ -oS) + \v , ; i \ 2 6(E„ + co)], (6) 



where the delta function 6(x) is taken as T/n{x 2 + T 2 ), with 
F = 0.01. The supercell technical is used to calculate the 
LDOS. 

In the following calculation, we take the hopping constant 
tjj to be unity for nearest neighbors and zero otherwise. The 
pairing potential V and the filling electron density n are cho- 
sen as V — 1.3 and n = 0.84 (hole-doped samples with doping 
6 - 0.16), respectively. The calculation is made on 48x48 lat- 
tice with periodic boundary condition and random distributed 
initial values of the order parameters are chosen. The 10 X 10 
supercell is used to calculate the LDOS. 

We summarize our main results in Fig. 1 , as seen, the H-T 
phase diagram is plotted. At zero temperatures two critical 
Zeeman fields hi - 0.14 and hi - 0.23, are revealed. The 
whole SC state is divided to be three regions, namely, uniform 



(f-wave SC state, ID FFLO state, and 2D FFLO state, respec- 
tively. The periodicity will decrease as the magnetic field in- 
creases in both ID and 2D FFLO states. As h > 0.28, the SC 
phase will be destroyed completely. The periodicity of the ID 
FFLO state also increases as the temperature increases, and 
it will transit to the uniform state as the temperature increases 
further. As a result, the range of the FFLO phase will decrease 
upon increasing the temperature. Near the SC transition tem- 
perature, only uniform cZ-wave phase and 2D FFLO phase was 
observed, with the transition field at about h = 0.225. 

The calculated order parameter amplitudes for various Zee- 
man fields h with the temperature T = 10~ 5 are shown in 
Figs.2(a)-2(e). As seen, for weaker magnetic field, the order 
parameter is uniform [Fig. 2(a)]. When the zeeman field in- 
creases, as we can see from Figs. 2(b) and 2(c), the SC order 
forms the stripe pattern. The order parameter is of nearly co- 
sine form with the periodicity of about 48 along x direction as 
h = 0.16. We have verified numerically that the periodicity is 
kept to be 48 for 0.14 < h < 0.175. And the periodicity re- 
duces to 24 as h increases (0.175 < h < 0.23). Here the finite 
size effect prevents us from obtaining solutions with period- 
icity not commensurate with the lattice size. As h increases 
further, the pattern changes to two-dimensional, with the pe- 
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FIG. 4: (Color online) The LDOS spectra for different phases. Panel 
(a) is the LDOS in the uniform phase with h = 0.13. Panels (b) 
and (c) are the spectra in the ID FFLO phase with h = 0.16 at the 
nodal line and at the site which the order parameter is maximum, re- 
spectively. The right panels are the spectra in the 2D FFLO phase 
with h = 0.235, where (d-f) are the spectra at the saddle point where 
two nodal lines intersect, the midsite between two neighboring sad- 
dle points, and the site where the order parameter has the maximum 
magnitude, respectively. The (blue) dotted line, (red) dashed line, 
and the (black) solid line are spin-up LDOS, spin-down LDOS, and 
whole LDOS, respectively. 




FIG. 3: (Color online) Plots of the FM order nij as a function of 
position for various Zeeman fields h. 



riodicity decreases as h increases, which can be seen clearly 
from Figs. 2(d) and 2(e). 

The spatial distributions of the FM order are shown in 
Figs.3(a)-3(e). The parameters are the same to those in Fig. 2. 
As seen in Fig. 3(a), in the uniform phase, the FM order is also 
uniform. Actually, the FM order competes with the SC or- 
der and is suppressed strongly by the SC order, as a result, the 
FM order could not survive at lower magnetic field (h < 0.12), 
and quite weak (« 0.003) as h — 0.13. We also checked nu- 
merically (not presented here) that the FM order will increase 
to about 0.05 in the normal state for the same magnetic field 
(h = 0.13). In the ID FFLO state, as seen in Figs. 3(b) and 
3(c), the FM order is largest along the nodal lines and is sup- 
pressed when the SC order parameter increases. The FM order 
reaches the minimum value as the SC order is maximum. The 
pattern also forms ID stripe but the periodicity is one-half of 
that of the order parameter. In the 2D FFLO state [Figs. 3(d) 
and 3(e)], the FM order forms the checkerboard pattern. Simi- 
lar to the case of ID FFLO state, the FM order is largest at the 
nodal lines and minimum as the SC order is maximum. The 
periodicity along the parallel direction is the same as that of 
the order parameter. While the periodicity along the diagonal 
direction is only one-half of that of the order parameter. 



We now turn to study the LDOS spectra. The LDOS [Eq. 
(6)] can be written as p,- = p,f +p,-^ . Here p,-f and p,| are respec- 
tively the spin-up and spin-down parts of the LDOS. These 
two parts are exactly the same if the Zeeman field is absent. 
In presence of the Zeeman field, the spin up LDOS shifts to 
left and the spin down LDOS shifts to right. In Figs.4(a)-4(f), 
we plot the two parts of LDOS separately to discuss the prop- 
erties of the LDOS. The whole LDOS spectra are also plotted 
so that the results can be compared with scanning tunneling 
microscopy (STM) experiments. 

The LDOS spectra in the uniform phase are shown in 
Fig.4(a). As seen, the spin-up LDOS shifts to the left with 
the mid-gap point locating at u> = —h. The SC coherent peaks 
shift to ±A(> - h. The spin-down LDOS shifts to the right 
with the SC coherent peaks at +Ao + h. Outside the gap we 
can seen the van Hove peak. As a result, the whole LDOS 
spectrum contains two stronger peaks at +(Ao + h) and two 
weaker peaks at ±(Ao-/i). The gap structure at low energies is 
"U"-shape. The density of states at zero energy p(0) increases 
linear with the external field, indicating that the quasiparticle 
excitations due to the magnetic fields. 

The LDOS spectra in the ID FFLO phase with h = 0.16 
are shown in Figs.4(b) and 4(c). Fig. 4(b) is for the site on the 
nodal line. We can see very sharp and strong peaks at the po- 
sition +h. The SC coherent peaks are suppressed and almost 
invisible. The peak at negative energy comes from the spin- 
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up LDOS, and the peak at positive energy are contributed by 
the spin-down LDOS. Taking into account the Zeeman shift, 
these in-gap peaks (bound states) at +h locate just at the mid- 
gap position. These bound states are due to the sign change 
of the order parameter across the nodal lines and are related 
to the Andreev reflections, similar to the mid-gap states in d- 
wave superconductors 11911 . The intensity of the in-gap peaks 
will decrease as the site moves away from the nodal line. As 
we can see from Fig. 4(c), at the site where the order parame- 
ter is maximum, the in-gap peaks are turned to be a hump at 
the mid-gap position for both spin-up and spin-down LDOS 
spectra. The SC coherent peaks are seen clearly. The mid-gap 
hump is so weak that it is concealed in the whole LDOS (pi) 
spectrum. We can see four peaks at the energies [+(Ao + h)]. 
And the spectrum of the whole LDOS is similar to that of the 
uniform phase while the gap structure at low energies is not 
"U"-shape but "V"-shape due to the presence of the mid-gap 
hump. 

At last we plot the LDOS spectra of the 2D FFLO phase 
in Figs.4(d)-4(f). Actually the features of the spin-up LDOS 
spectra are studied intensively in Ref. [15]. There are two 
kinds of Andreev bound states. One is due to the sign change 
of the order parameter across the nodal lines. The second is 
essentially localized at the saddle points. The order parame- 
ter is suppressed strongly in a intersecting region, which pro- 
duces a potential well for a quasiparticle and thus generates 
two finite-energy andreev bound states. As a result, at the 
saddle points, four in-gap peaks exist in the spin-up LDOS 
spectra [Fig.4(d)] at the energies -0.575, -0.385, -0.08 and 



0.105. And mid-gap peaks exist between two neighboring 
saddle points [Fig. 4(e)]. At the site where the order param- 
eter is maximum, the LDOS spectrum [Fig. 4(f)] is similar to 
that of the ID FFLO state [Fig. 4(b)] and that of the uni- 
form phase, namely, if the van hove peaks and the weak peak 
caused by the mid-gap hump are excluded, there are only four 
peaks left, locating at +(Ao + h), contributed by the spin-up 
and spin-down LDOS, respectively. 

We have shown the LDOS spectra of the three different 
phases. As seen in Figs.4(a)-4(f), the spectra are quite differ- 
ent and the spectra in each phase have their distinctive features 
as we discussed above. Thus they can be easily detected by 
the STM experiments and can be used as signatures to probe 
the FFLO states. 

In summary, based on a BCS-type model and BdG equa- 
tions, we studied the phase transition induced by the external 
magnetic field. The phase diagram is mapped out and the tran- 
sitions from uniform phase to ID FFLO state, and ID FFLO 
state to 2D FFLO state are revealed. We also calculate the 
LDOS to discuss the signatures of the three phases, namely, 
the LDOS spectra in the uniform phase will contain four peaks 
due to the Zeeman shift. In the ID FFLO state, the LDOS 
spectra show mid-gap states due to the Andreev reflection. In 
the 2D FFLO states, four in-gap peaks are revealed at the sad- 
dle point due to two kinds of Andreev bound states. 
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